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The effective thermal conductivity of composites with ellipsoidal fillers is analyzed by using a homogenization 


method that is able to represent the microstructure precisely. In this study, various parameters such as the volume 


fraction, shape, and distribution of the filler are quantitatively estimated to understand the mechanisms of heat 


transfer in the composite. First, thermal boundary resistance between resin and filler is important for obtaining 


composites with higher thermal conductivity. Second, the anisotropy of the effective thermal conductivity arises 


from contact between filler in the case of ellipsoidal filler and produces lower thermal resistance. Finally, the fil- 


ler network and thermal resistance are essential for the heat transfer in composites because the path of thermal 


conduction is improved by contact between neighboring filler particles. 
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Introduction 


The thermal conductivity of resin is typically around 
0.15 to 0.3 W/(m-K), which is much lower than that of 
metals and ceramics. Resin has thus often been used as 
an adiabatic material and for reinforcement due to its 
advantages of lower heat conduction and higher strength. 
In the last two decades, demand has grown for resins 
with higher conductivity and better heat dissipation cha- 
racteristics in the field of electronics because of the better 
formability of resins compared with metals and ceramics 
[1, 2]. However, there are limits on how much the con- 
ductivity can be increased. To obtain higher conductivity, 
modification of composites has been performed by the 
addition of filler, which is a metal and ceramic of higher 
thermal conductivity. However, the mechanism of heat 
transfer in composites is quite complex because of the 
following factors related to obtaining effective thermal 
conductivity (ETC) [3-5]. 

1. Thermal conductivity of resin and filler 

2. Volume fraction of dispersed filler material 

3. Shape and size of filler 

4. Dispersion condition of filler particles 

5. Thermal boundary resistance between resin and fil- 


ler particles [6-8] 
6. Temperature distribution around filler particles 
7. Orientation of filler particles for thermal conductivity 


There are many methods for analyzing the ETC of 
composites [6-11]. However, it is quite difficult for con- 
ventional analytical models to model the above-men- 
tioned factors perfectly. Direct numerical simulation, 
which can represent the microstructure of a composite 
with fillers of various shapes and volume fraction is re- 
quired for precise thermal analysis. In our previous stu- 
dies, we developed a homogenization method for thermal 
analysis of porous media [12] and packed beds [13] with 
complex structures. The method employs a three-dimen- 
sional finite element method that enables calculation of 
the ETC, making this approach very useful because it is 
able to consider changes in microstructure precisely. For 
example, the geometric anisotropy of the microstructure 
depending on pore shape and configuration was ex- 
amined precisely, and the importance of microstructure 
was shown quantitatively for heat transfer in porous me- 
dia and packed beds. Accordingly, this homogenization 
method is attractive for obtaining more quantitative in- 
formation about the above-mentioned factors in compo- 


sites [14, 15]. However, little research has focused on 
perfect modeling of the volume fraction, shape, and dis- 
tribution of the filler in composites. In particular, expe- 
rimental estimation of thermal boundary resistance be- 
tween resin and filler particles has been difficult although 


it is caused by incomplete interfacial contact such as void. 


In this study, we investigate thermal behavior in compo- 
sites while simultaneously considering factors 1, 2, 3, 4 
and 5 described above by means of our homogenization 
method. The resistance between fillers might be impor- 
tant. Moreover, although this method can also simulta- 
neously consider factors 6 and 7 and different size of 
fillers, this study focuses on factors 1 to 5 which are re- 
lated to heat transfer in composites. Finally, the ETC of 
composite is compared to that of porous medium with 
ellipsoidal pores. 


Model 


Homogenization method 

To analyze the composites shown in Fig. 1(a), the 
simple periodic structure shown in Fig. 1(b) is assumed. 
The unit cell of the periodic structure consists of two 
domains: solid 1 (resin), Q,), and solid 2 (filler), Qu, as 
shown in Fig. l(c). Z” denotes the interface between the 
two domains. 


(a) Composites (b) Periodic 


(c) Unit cell 


Fig. 1 Schematic diagram of homogenization method 


The microscopic scale (filler size, /) in the periodic 
domain Q is small compared with the characteristic 
length L at the macroscopic scale. 

The multiscale periodic heat conduction problem un- 
der steady-state conditions without heat generation can 
be mathematically expressed as 
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where / and T are the thermal conductivity and tem- 
perature field, respectively, and h is the coefficient of 
thermal boundary resistance between the two domains. It 
is caused by the incomplete contact. The Biot number is 


given by 
Bi=—. (4) 


sl 

Finally, the ETC is obtained as a homogenized para- 
meter by multiscale asymptotic expansion. Details of the 
calculations are described in the packed bed analysis of 
our previous study [13], with gaseous and solid phases in 
the reference corresponding to solids 1 and 2 in the com- 
posite, respectively. 

In the present study, finite element meshes of simple 
cubic (SC), body-centered cubic (BCC), and face- 
centered cubic (FCC) configurations are created as filler 
positions in the composites by using a tetragonal mesh by 
the Delaunay triangulation algorithm [16]. Ellipsoidal 
filler particles with b = c > a are distributed in the cell as 
shown in Fig. 2. ETC is calculated for composites with 
ellipsoidal fillers of a/c = 0.25, 0.50, 0.75, and 1.0 for 
various Biot numbers and volume fractions p, and the 
mechanisms behind the manifestation of anisotropy of Ax 
and Ay (= A,) are examined. 
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(b) View of composite 
with ellipsoidal filler 


(a) Size parameter 


Fig. 2 Unit cell of composite containing ellipsoidal filler 


Materials 

The thermal conductivity of the epoxy resin, Ası is 
kept constant (A,; = 0.2) in this study. The thermal con- 
ductivity of the filler, Ay, is set to 5, 30, and 170 
W/(m-K), which corresponds to the conductivity of metal 
and ceramic (AI,O; and AIN). 


Results and Discussion 


Effect of filler fraction and filler shape on ETC of 
composite 

Figures 3, 4, and 5 show the effect of filler shape with 
SC geometry on ETC for the thermal boundary resistance 
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when the volume fraction is p= 0.1, 0.3 and 0.5 and 
the conductivity of the filler is A = 5, 30, and 170 
W/(m:K). Solid symbols and open symbols indicate val- 
ues for A, and 4y, respectively. When the filler particles 
are isolated or dispersed well at lower faction, ETC 
strongly depend on the thermal resistance between the 
two domains. For example, if the thermal resistance is 
higher (lower Biot number, Bi < 0.1), thermal conduction 
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Fig. 3 Effect of thermal resistance and volume fraction on 
ETC for the SC model at A, = 5 
through the interface between resin and isolated filler can 
be ignored. If the heat transfer behaves like a porous ma- 
terial, then the composite acts as an adiabatic material. 
As aresult, ETC is similar to the conductivity of the resin 
because the resin is a continuous phase. However, if a 
medium Biot number is used, which means moderate 
thermal resistance, then the thermal conduction through 
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the interface between the two domains is improved. The 
ETC gradually increases with Biot number, and then le- 
vels off. Because thermal conduction through the inter- 
face becomes perfect in the case of higher Biot number 
(Bi > 10000), the saturated ETC is proportional to the 
volume fraction and the conductivity of filler. Therefore, 
the conductivity of filler becomes more important for 
improving the thermal conduction of the composite as the 
Biot number becomes higher. 
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Fig. 4 Effect of thermal resistance and volume fraction on 
ETC for the SC model at A, = 30 

When the filler shape is flatter, the anisotropy of the 
ETC becomes more remarkable. The anisotropy is not 
particularly high because filler particles are not in contact 
with neighboring filler particles at p= 0.1. However, as 
the fraction becomes larger, thermal conduction improves 
in only the y-direction in the case of flatter filler. When 
the fraction reaches 0.5, filler particles are in contact in 


y-direction even with nearly spherical filler (a/c = 0.75). 
Accordingly, the anisotropy of the ETC depends on the 
volume fraction and oblateness of the filler particles. 
However, the thermal conduction from the pathway in 
the x-direction that passes through the solid 1 (resin) 
phase decreases. For the same fraction, it is easier for 
flatter filler to produce closed pathways. For example, in 
the case of SC geometry with p= 0.3 and a/c = 0.25, the 
pathway in the x-direction through the solid 1 phase is 
almost blocked by the solid 2 (filler) phase. In the case of 
lower Biot numbers (higher thermal boundary resistance), 
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Fig. 5 Effect of thermal resistance and volume fraction on 
ETC for the SC model at A = 170 
interruption of the pathway by solid 2 acts adiabatically 
in the x-direction although the ETC in the y-direction 
increases slightly due to contact between filler particles. 
Finally, the ETC of composites with flatter filler exhibits 
anisotropy between the x- and y-directions when the Biot 
number is lower than 10. The anisotropy at lower Biot 
number strongly depends on the filler shape. However, in 


the case of higher Biot numbers, the thermal resistance 
through the interface is improved in the x-direction, and 
the anisotropy almost disappears due to the perfect ther- 
mal conduction. When filler particles are in contact with 
neighboring particles, the anisotropy increases. This 
means that the connection network of the filler is impor- 
tant for obtaining higher thermal conduction. In particular, 
flatter filler exhibits stronger anisotropy at the same vo- 
lume fraction. ETC in the y-direction for flatter filler at a 
larger fraction becomes almost constant regardless of 
Biot number. Thermal resistance does not affect ETC in 
y-direction because most of the interface of flatter filler is 
nearly parallel to the y-direction. Moreover, it becomes 
easier to establish a filler network as the volume fraction 
increases. However, composites with nearly spherical 
filler do not exhibit anisotropy for every Biot number. 

A similar tendency is obtained for filler with higher 
conductivity, although the absolute values are higher. The 
filler shape, thermal conductivity, volume fraction, and 
thermal resistance through the interface are therefore 
essential for correctly estimating ETC. Moreover, if ETC 
is measured experimentally, quantitative prediction of 
thermal resistance, which is caused by incomplete con- 
tact between resin and filler particles, might become 
possible. 

Effect of filler configuration on composite ETC 

The effect of filler distribution on ETC is analyzed 
because thermal boundary resistance is important for 
actual composites with flatter fillers. ETC for different 
filler fraction and filler shapes is calculated for three con- 
figurations (SC, BCC, and FCC) at A, = 30, and at Biot 
number Bi = 1, 100, and 10000 (low, medium, and per- 
fect thermal conduction through the interface) as shown 
in Figs. 6, 7, and 8. In these figures, solid symbols and 
open symbols indicate results for A, and ày, respectively. 
For smaller a/c, the anisotropy of ETC becomes larger at 
around p= 0.2-0.7 for each geometry when the Biot 
number is low (1) and medium (100). For example, in the 
case of low volume fraction and low Biot number, ETC 
is smaller. When the flatter fillers come into contact with 
each other, ETC in the y-direction suddenly increases at 
the volume fraction where a pathway of filler with higher 
conductivity penetrates through the continuous resin 
phase. However, ETC in the x-direction remains low in 
the case of low Biot number and flat filler. An increase in 
the ETC in the x-direction appears suddenly upon pene- 
tration in the x-direction. Two sudden increases in ETC 
thus appear in the case of flatter filler. The difference 
between the two depends on the oblateness of the filler. 
Moreover, the slope until the penetration points depends 
on the Biot number. For example, the slopes are quite 
low and gentle for Bi = 1 and Bi = 100. Finally, the two 
discontinuities caused by the penetration disappear in the 
case of perfect thermal conduction. 
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The effect of filler configuration on ETC is larger 
when the volume fraction is around 0.2 to 0.5. Because 
the higher conductivity solid 2 phase penetrates the lower 
conductivity solid 1 phase, pathways formed by a net- 
work of solid 2 phase are essential for increasing the 
thermal conduction in the composite. For example, ellip- 
soidal fillers can make contact with neighboring filler 
particles easily even for the same volume fraction. Con- 
tact between filler particles in the FCC and BCC geome- 
tries is slightly more difficult than in the SC geometry 
because the distances between the filler particles are dif- 
ferent. The difference in connection level caused by flat- 
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Fig. 6 Effect of volume fraction on ETC for the SC model at 
Ag = 30 

ter fillers affects the anisotropy as well because connec- 
tion becomes more difficult for flatter and horizontal 
fillers in the vertical direction than in the horizontal di- 
rection. 

However, in the case of higher Biot number, a linear 
relation between the volume fraction and ETC is ob- 
tained for each geometry and filler shape. Because ther- 
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mal conduction through the interface becomes perfect, 
ETC becomes proportional to the volume fraction and 
conductivity of the two components. Finally, the filler 
configuration and shape must be taken into account for 
thermal analysis of composites when the thermal resis- 
tance is not perfect. 
Comparison with ETC of porous medium with ellip- 
soidal pores 

In previous study [12], ETC of porous medium with 
ellipsoidal pores was analyzed and the effects of parame- 
ters such as porosity, pore shape, pore distribution, and 
temperature of the porous medium on the conductivity 
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Fig. 7 Effect of volume fraction on ETC for the BCC model 
at Aso = 30 
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Fig. 8 Effect of volume fraction on ETC for the FCC model at 
As2 = 30 


were studied to clarify the mechanisms in complex pore 
structures. In the case of composite with fillers, thermal 
conductivity of dispersed phase (filler) is higher than that 
of continuous phase (resin). Accordingly, ETC of com- 
posite is improved by thermal conduction through the 
pathway of dispersed phase. As mentioned above, filler 
network, which is a connection between neighboring 
filler is influential to the ETCs of composite. 

On the other hand, in the case of porous medium, 
thermal conductivity of dispersed phase (gas) is lower 
than that of continuous phase (solid). Even if ellipsoidal 
pore has contacts with neighboring pores, the pathway of 
continuous phase of higher conductivity is unchanged. In 
this case, the ETCs are relatively constant. However, 
when the pathway for thermal conduction is completely 
blocked by gaseous phase of lower conductivity, the 
ETCs drastically decrease. Accordingly, cross-sectional 
pathway of the solid phase plays more roles than the pore 
connections. The effect of thermal radiation on the ETCs 


could be estimated for porous medium, meanwhile ther- 
mal radiation effect on the composite is insignificant. It 
was observed that thermal radiation is dominant at high 
temperatures and for large pore sizes, and the anisotropy 
introduced by pore shape could be neglected. Based on 
the aforementioned, the ETCs of materials with ellip- 
soidal pores or ellipsoidal fillers showed different beha- 
viors under homogenization method which indicate the 
microstructure and the conductivity of each phases. 


Conclusions 


A homogenization method that is able to consider heat 
transfer with thermal resistance was applied to thermal 
analysis of composites with ellipsoidal filler. The effect 
of microstructural aspects such as filler configuration, 
shape, and volume fraction on ETC was calculated. First, 
ETC depends on Biot number, which is related to thermal 
boundary resistance, with the value drastically increasing 
at Bi = 10 to 1000. In the case of lower Biot number, the 
composite behaves like a porous material. Second, larger 
anisotropy of ETC arises from flatter fillers and higher 
thermal resistance because pathways through the differ- 
ent phase with lower conductivity are blocked. Thermal 
resistance and contact between filler particles are impor- 
tant for heat transfer in the composites. Moreover, quan- 
titative data considered for contact between the filler par- 
ticles show that the anisotropy of ETC occurs easily 
when ellipsoidal filler is placed in the SC configuration 
compared to FCC and BCC configurations. Finally, the 
proposed three-dimensional method is a useful tool for 
precisely modeling heat transfer with thermal resistance 
in composites because it can estimate the microstructure 
of filler shape and configuration. 
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